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We have previously applied the complete first-order adiabatic 
correlation function of Dalgarno and Lynn to the problem of low-energy 
elastic positron scattering by hydrogen and helium atoms. This approach 
is now extended to yield rigorous lower bounds on the scattering phase 
shifts in the case of hydrogen and "quasi-bounds" for helium. In addition, 
the positron annihilation rate in helium is re-evaluated, and lower values 
are now found, with the enhancement factor over the Dirac rate varying 
between 1.5 and 1.8. Some comparisons are made between the present method 


and other recent work. 
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In two recent papers 1 * 2 , the problem of low-energy elastic 
scattering of positrons from simple atoms (hydrogen and helium) has 
been treated in a modified adiabatic approximation. The method con- 
sisted in assuming that the optical potential 3 is well represented by 
the position-dependent second-order energy shift in the ground state 
of the target atom, due to the electric field of the positron 4 . Two 
novel features were introduced in these calculations. Firstly, the 
exact second-order potential 5 was used in place of an approximate form 4 . 
Secondly, suppression of the short-range parts of the potential was 
accomplished by the use of a semi-empirical energy independent parameter, 
which served to reduce the spherical (monopole) distortion of the atom 
which would otherwise dominate as the positron approached the nucleus. 

In the case of pos itron -hydrogen scattering 1 , 90$ suppression of 
the monopole was found to give essentially exact agreement with the 
definitive S-wave results of Schwartz 6 from zero energy up to the posi- 
troniun formation threshold. Higher partial waves gave smaller phase-shifts 
than the best present estimates 7 , presumably due to the inclusion of 
insufficient amounts of virtual positronium in the wave function. 

The pos itron -helium calculation 2 was repeated for various values of 
monopole suppression, and suffered frctn the use of the shielded hydrogenic 
approximation for the atomic ground-state. But since no exact results 
exist its correctness cannot yet be evaluated 8 . The only existing experi- 
ment 9 measures the momentum transfer cross-section near the positronium 
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formation threhold, and seems to disagree strongly. 

The results for hydrogen 1 seem to imply that the adiabatic wave 
function contains, to a considerable extent, the significant physics of 
low-energy positron scattering. Nevertheless, it is disturbing that a 
semi-empirical parameter is required, and also that no extremum principle 
or bound is contained in the results. The purpose of the present work 10 
is to use the adiabatic wave function in such a way as to obtain rigorous 
lower bounds to the phase shifts for e + - H scattering and "quasi-rigorous" 
lower bounds in the e + - He case. 

In Section II, the lower bound principle of Gailitis 11 is adapted 
to the adiabatic type of trial function for e + - H scattering. In 
Section III, a simple scale change gives the corresponding equations for 
e + - He scattering, and the e + annihilation rate in helium is also calcu- 
lated. Results for these two cases are given and discussed in Section IV. 
In Section V we clarify the relations between the present method and other 
recent work. Two Appendices contain details of the calculation. 

II. Formulation and Positron -Hydrogen Scattering. 

We are interested here in obtaining approximate solutions to the 
positron-atom elastic scattering problem. The energy will be restricted 
so as to forbid any re-arrangement or inelastic processes, the lowest of 
which is real positronium formation, occurring above 6.8eV in hydrogen 
and 17.8eV in helium. To apply the lower bound principle of Gailitis 11 , 
it is sufficient in this energy range to use a trial scattering function 


of the form 
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Here, x is the position vector of the positron, r represents the atomic 

<^v_/ 

electrons, 0 is the ground-state wave function of the target atom and X 
is the scattering function for the positron of momentum k, which must 
approach the correct asymptotic form as x -* <». The distortion or polar- 
ization of the target by the incoming positron - the closed-channel part 
of the wave function - is approximately represented by §, which satisfies 
the conditions 12 


<p(4) $ = O 


for all x. 




(2b) 


Then the two functions x and § are to be determined. The resulting 


form of x yields partial -wave phase shifts § (k) which satisfy the 

J_i 

inequality tan 6 L (k) g tan 6 L (k) (exact), for all energies below any 
resonances that might occur. 

The Hamiltonian for e + - H scattering is 

h = h* + h*. + y 

where (- 1 ^ = - [ + 2 A 1 > 


and 


VC 
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where we have set h 2 / 2m e = 1, e 2 - 2, with energies in Rydbergs and 
lengths in units of a , the Bohr radius. The correlation function is 
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<£> (A S ) = F0&) &&>%'> 


where F is to contain short-range non-adiabatic modifications to the 
adiabatic correlation function G which satisfies the first-order equation 5 

[ & , HJ 0^ = [V- <y>l 0^- (5) 


[Throughout this paper the bracket notation will represent expectation 
value in the ground state, ie 


<y> = (6) 


From Eq. ( 5 ), the function G is determined to within an arbitrary additive 
function of x, which may be chosen to make <G> vanish. With this choice, 
Eq. ( 4 ) is consistent with Eq. (2a). 

If we use standard variational methods 13 anri allow X an( i E to undergo 


r 
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free, independent variation, we obtain the following two coupled 
equations : 




(7a) 


<^Gr [ H-El [7- +■ FG-l') = O; 


(7b) 


where E _ E q + k 2 and (H^ - E Q )0(r) = 0. Simple operations, including 

commutation of H and G and use of Eq. (5), yield the following coupled 
r 

differential equations for x and F: 


[vm 2 - va* = f > 

[w (v z tl0 + v 2 - V 3 - w - ] F = \4 % . 


The "potentials" appearing in Eq. (8) are defined as follows: 
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V t (*) = <& P 'v> , (r-' A3 )i 


W (*") - 
h / (*0 = 
V„ (3C) - 


- <fr\‘&>J 

<G- 2 >; 

d W 
doc “ 


(9) 


Using the exact analytical form of G, given in elliptical coordinates 
by Dalgarno and Lynn 5 , we have numerically evaluated 14 all the functions 
defined in Eq. (9). The evaluation is outlined in Appendix A, and in 
Table I the potentials are tabulated. 

The limiting forms of these potentials can be found analytically for 
small and large x, and are useful for starting the solutions of Eq. (8) 
and also to check their numerical evaluation. (See Appendix A) We obtain, 
for x -► o 


V t « z Ik-'] 

V z ^ -1 -i- o o») 


N 

w ^ 




fc [1-2*1 , 


(10) 


V * 
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and for x CO 

Y 1 and W -* 0 exponentially. 


V 2 ^ — [9/(2x 4 ) + * ] 

V 3 ^ - [2i?> /( 2X ' 7 ) +■ , ' 7 ^/ acq * 3 

h/ /W [43/(8X-») 4- |07/(8X 4 )f ] 

^ [4$/C** 9 ) ]. 


(ii) 


The solution of the coupled Eqs. (8) is outlined in Appendix B, and 
the results are discussed in Section IV. In Fig. 1 we have plotted the 
S-wave zero-energy solutions for the two functions X and F - x, as well 
as the "suppression factor" T(x) - F/x- One can see the natural occurrence 
of the expected short-range suppression of G. An unexpected small enhance- 
ment also occurs for x > 4. 

III. Positron-Helium Scattering 

To extend the results of the preceding Section to the case of positron 
scattering from the ground state of helium, we interpret the correlation 
function of Eq,. (4) as depending on two electronic coordinates, ie 


$( 


Jl< . A-i 


S) _ FG5) G-Ov-H*) 


( 12 ) 


where 0 is the properly antisymmetric ground-state atonic wave function. 
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and G is symmetric in the electronic coordinates. 
The Hamiltonian for e + - He scattering is 


H = H 12 + * V, 


(13) 


where 


h„ - -h** V+ h * k\ 


A: 


Hor - 


- V, 


and 


V = V a) v v‘°. 


where 


, <-> „ t 

V - 2 [ 




1. 


The assumption that G is the first-order adiabatic correlation function 
requires 2 


[ = [V- <y>] $ fa J A 2 ') } 


( 14 ) 
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where the bracket notation of Eq. (6) is extended to include integration 
over r^ and r 2 . 

The formal equations to be solved are identical to Eq. (j) where 
E _ E q + k 2 and (H l2 - E Q )0(r 1 , r 2 ) _ 0. Using Eq. (l4) one obtains the 
following coupled equations: 

~ <G-V> F, (ua) 


(15b) 

+ ?*<&*>• Y*If = <&v>7^ 


If it were possible to solve Eq. (l4) exactly, and to evaluate all 
the bracketed expressions in Eq. (15), the resulting phase shifts would 
satisfy the rigorous lower bound theorem. Since no exact solution for 
the helium atom ground state is known, one must use an approximate form 
for 0(r 1 , r 2 ). We have used the simple, shielded hydrogenic form 0(r 1 , r 2 ) 
T l r ( ;r i ) \Jf(r 2 )S, where S is a singlet, antisymmetric spin function, and 
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(p M - 



(16) 


where g is an effective charge. Since this function is not an exact 
eigenfunction of H 12 , the procedure leading to Eq. (15) is itself not 
really consistent. To the extent that these inconsistencies may be ignored 
however, the solutions of Eq. (15) will yield lower bounds to tan we 
call these "quasi -bounds". 

Using this form of <p, it is consistent with Eq. (lh) to assume 
independent distortion of the two electrons. 


C-( 


A t 




^0 - Or + G" + G~z- (1T) 


Eq. ( 1 ^ 4 - ) then simplifies to the form 


H ^ y 2 ] ■= [ <v w > , v (i) J , 

(z - u i) 


(18) 
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■where bracketed quantities now refer to single integrals, ie 


<v a ’> ^ fa </>&) V W O' M . 


( 19 ) 


Again, one must set <£h> = 0 to satisfy Eq. (2a). With this condition 
satisfied, one can immediately re-write Eq. ( 15 ) in terms of single 
integral quantities : 


- 2V,F, (“a) 

+ V 2 -Vj-W-V Wd ^] F -- V/A. 

(20b) 


All the "potentials" used in Eq. (20) have the formal definitions of 
Eq. ( 9 ), but with $(r) replaced by \|r(r). The factors of 2 appearing in 
Eq. (20a) are due to the superimposed effects of the two electrons. In 
Eq. (20b), there is an additional term, -NV^ which did not appear for 
hydrogen, Eq. ( 8 b). This comes from the third-order term in Eq. ( 15 b), 
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<c-“v> = <(&,*&,') (v"v")> 

= <5,‘v"> * <Ci‘v“’) 

* <<?,’> <y‘"> • <s’7< v "> 

+ Z <&1> < & i vW ) 4-2 <& a ')<G 1 V u )_ 


The last two terms vanish, and the net result is 

<frV> = 2 [V 5 + NV,]_ ( 

It has been noted 2 that a simple change of scale converts Eq. (l8) 
into its hydrogen equivalent, Eq. (5}. Specifically, let y - pr and 

ro /v* 

v _ gx. Then Eq. (l8) becomes 
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One sees immediately from Eqs. (23) and (5) that 

(3 Gfe - G h CIO- 


(24) 


Also, since V Hg (r, x ) = pVgCy, v), and ^(rJcPr = 0^(y)d 3 y one can relate 
all the helium potentials of Eq.. (20) to the hydrogen potentials of Eq_. (8) 
as follows : 


K / (pc) — ^ p f ^ ( u X 

V„ (*) — > f' V " ^ 
v (oo — - > w (o . 


(25) 
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The simplest way to make use of these scaling laws is to change variables 
throughout Eq. ( 20 ): 


^ I'V - p- 


(26a) 


[* W ( V 1 +Ji! - (3% ^ ^dxr 1 


(26b) 


- y; 'V ■ 


Here k 1 _ k/g, V 2 = V 2 , and all the potentials are the hydrogen functions 
±1 terms of v. These equations can be treated just like those appearing 
in the hydrogen problem, as discussed in Appendix B, and the results are 
discussed in Section IV. As indicated previously 2 , the value 3 = Jg 
which minimizes the energy expectation value gives an unsatisfactory value 
(l.ll) for the polarizability of helium. This quantity appears in the 
present theory as lim [ 2 x 4 V 2 (x)] = 93 4 * With only a slight change in 

X -* OD • 

the energy, one can adjust the polarizability to the experimental value 2 
of 1 . 376 , by setting g - I.5992. This value has been used throughout the 
present work. 

We have previously 15 used the non -variational e + - He wave functions 


to calculate Z^, the number of atomic electrons which are effective in 
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annihilating a positron of momentum k. (For planewave positrons 
and omitting atomic distortion, _ 2) . In terms of our type of 
wave function, 

i-z K = 

_ Q»3C <|/ 2 £» [ I y i- F G- I Fl *<&*>] 


(27) 


Following Ref. 15, we assume that only the L _ 0 scattering function 
diffeis significantly from its planewave form, and that for L > 0 the 
ratio of closed-channel to open-channel functions is well approximated 
by the zero-energy suppression function T. With these approximations, 
one can write 




% G- + 9c 2 



* 4 


Ur e' Zi ''V MM[i-4^(p 


(28) 


c 



* V 
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where 


M<y) - f i +• ~t &(vr.ir)]^ <■ T 


WO*') 


Here v _ gx, U q and g Q are the S-wave parts of X and F (see Eq. Bl). 

The correlation functions used here [G(v, v) and W(v)] are scaled 
appropriately for helium. At k = 0 the second integral in Eq. (28) 
vanishes since no higher L-values contribute. For k + 0 the assumptions 
made here should be quite good; in Ref. 15 we have estimated that an 
error of less than 2 <f is incurred in the values of Z, . In Table II the 
results are given for values of the energy up to the positroniun thres- 
hold, in comparison with the larger values of previously obtained 15 . 

The difference between the present results and those obtained 
previously is traceable partly to the fact that the present calculation 
gives less attraction than did the non -variational method 15 , and hence 
the positron is pulled inward to a smaller degree. A larger effect is 
the difference between the "correlation function" M defined in Eq. (28) 
and the monopole - suppressed function Y used in Ref. (15)* These two 
functions are compared in Table 113^ and one can see that, except at very 
small distances from the nucleus, M < Y. To examine these two effects 
quantitatively, we have re-computed ^ at zero energy using both M and 
Y, as well as both variational and non -variational scattering functions. 

By far the largest effect is produced by replacing M by Y, outweighing 
the effect of interchanging the two forms of U by a factor of 3- Since 
neither M nor Y is exact, this result constitutes a warning against taking 
the annihilation results too seriously. It is interesting, however, that 
the present calculation agrees fairly well at zero energy with the expert- 
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mental result of Falk et al ls for thermal positrons of Z _ 3.78 ± O.17. 

k 

IV. Scattering Results and Discussion 

The e + - H S-wave scattering phase shifts 6 q were computed numerically 
as discussed in Appendix B, where the method of treating the asymptotic 
behavior of the closed- channel is discussed. This involves inward 
integration from a point x q in the asymptotic region to an intermediate 
point x. In Table IV the sensitivity of tan 6 to the choice of x q 
and x is shown for low and high-energy positrons, confirming the stability 
of the numerical solutions . 

In Table V the e + - H S-wave phase shifts (and, at k _ 0, the scattering 
length) are shown, along with the results of Schwartz 6 . To assess the 
quality of the present work, we compute the quantity the difference 
between tan 6 q and the value of tan obtained with an undistorted atomic 
wave function. This "Hartree" result is also given in Table V. The 
figure of merit Q(k) ..= A (variational ) /a (Schwartz), measures how much of 
the distortion or polarization has been accounted for by our method. It 
is most interesting that Q « 90f o , almost independent of energy. The 
lower bound results of Hahn and Spruch 17 have somewhat lower Q-values, 
as shown in the Table. We consider it a success of the present method 
that the results compare well with the many-parameter, many- l variational 
results of Ref. 17- Their method has the considerable advantage, however, 
of being capable of systematic improvement, by the addition of more terms. 

The results can also be convincingly extrapolated, and give good agreement 
with those of Schwartz 6 . 

In Table VI the S-wave scattering results for e + - He are presented. 
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Along with the "quasi -rigorous" lower hounds obtained here, we show 
the "Hartree" results, obtained fran a two-term analytic form 2 
approximating the helium ground-state potential, and the non -variational 
results 2 . Since no definitive calculation exists in this case, we have 
computed the figure of merit Q assuming that 6^ is exact. These 
Q-values are smaller than those obtained for hydrogen, but are also 
nearly constant. It is not possible to make any further statements about 
the "correct" values for the S-wave e + - He phase shifts from the com- 
parison of the quasi -variational and non-variational results. 

The possibility that a bound e + - e - P exists has been eliminated 
by an extensive variational calculation 18 which showed that for a 
"positron" of mass m > 2. 625 such a bound system does occur. As a 
further test of the present method, we have calculated this critical mass 
by looking for the appearance of a node in the S-wave zero-energy scattering 
function, corresponding to a scattering lenght a = ± oo. One modifies 
Eq. (8) by dividing all terms involving V 2 by the "positron" mass m. This 
is equivalent to multiplying the potentials V 1 , V 2 , V 3 by m, and adjust- 
ing the asymptotic conditions in an obvious way. We find the critical 
mass to be 2-79> and since our method gives an upper bound this agrees 
with the result of Ref. 18. 

The present method is not expected to be as good for L > 0 as for 
L - 0, since even the non-variational results 1 underestimate the P- and 
D-wave phase-shifts at the higher energies. This defect seems to be 
associated with a lack of sufficient virtual positronium in the wave 
function. Support for this is the success of Bransden' s P-wave results 19 
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at high energy. More recent results 20 indicate that just above the 
threshold for positronium formation, P -waves dominate over S-waves in 
the pick-up process, implying that below threshold neglect of positronium 
is less serious for S-waves. In any case, we show our L _ 1, e + - H 
results in Table VII, where they are compared with the lower-bound and 
extrapolated results of Kleinman et al 21 . 

V. Comparison with Other Methods 

We will now describe briefly some other methods which resemble the 
present work, hoping that some clarification will result. The field of 
low-energy positron scattering suffers from some lack of precision in 
terminology and some overlapping of approximations, which may make it 
difficult to evaluate the work reported here without a description of 
the competing methods. The following short review is not complete, but 
it may aid the orientation of the reader. 

1. The non-variational (polarized orbital) method 1 uses the form 
of scattering function given in Eqs. (l) and (4) with the 
additional restriction F(x) = x(x) . Since it is non-variational, 
however, the single unknown function x is determined by simply 
projecting the SchrCdinger equation onto the open-channel 
part of 


<( [H-Elj>&]> % — O. 


(29) 
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This leads to the differential equation 

[v 2 +f\% - lv,*v,]?Uo : 


whose solutions gives the phase shifts. In the absence of 
an applicable lower bound principle it is difficult to gauge 
the error in the phase shifts obtained from Eq. (30). 

2. The same polarized orbital wave function can be employed 
variationally . The same procedure which leads to Eq. (7) 
when X and F are independently varied, leads now to 


L 1 V- - 0 


(31) 


since F = x> This gives the differential equation 




Vj fVj <-W 




(32) 
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whose solutions provide approximate phase shifts which 
are lower bounds to the exact ones. Unfortunately, these 
approximate phase shifts are very low (a — - 1 . 0 ), since 
the large, short-range repulsion represented by W, and 
the reduction in V 2 caused by the factor (l + W )" 1 more 
than overcome the attraction added by V 3 . It is likely 
that removing some or all of the monopole part of G, as 
was done non-variationally in Ref. 1, would improve the 
results considerable. If G were replaced by G - (l - a )G , 
the potentials appearing in Eq. ( 32 ) would be modified as 
follows : 

N — ■> is/ - (i-^Wo 
w — > w — (1-^^° 









— > N/^ — (1 ^ V 2C , 

v 3 — * v 3 - (a - <>0 -2 <j-*)p 


(33) 


where, for example, N = <G 2 >, and P = <G (G - G ) (V-V )>. 

00 o o o 

In Eq. (33 ) > "the parameter a represents the fraction of 

adiabatic monopole distortion being retained; a - 0 is 

what we have called 1 "full monopole suppression". The best 

value of 01 , for each L and k, is to be determined variationally. 
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Eq. (32) is integrated for several values of a, and the 
value which maximizes the phase shifts is retained. The 
monopole terms (with subscript 0) are easily evaluated 
in spherical coordinates, while the remaining terms are 
kncwn in elliptical coordinates, except for the last 
term in V 3 . The quantity P is difficult to compute in 
either set of co-ordinates, but it might be approximated 
using the multipole expansions for G and V. This method 
has not been used up to the present time. 

3. Callaway et al 22 have made some modifications to Eq. (32) 
and arrive at another optical potential which they call the 
extended polarization potential. Although the results 
obtained are poor, the new potential has considerable 
physical interest. 

To derive the result of Ref. (22), one first neglects 
all terms which are formally of higher than second order in 
the perturbation V. This eliminates from Eq. (32) all re- 
ference to N and V 3 , and leaves us with Eq. (30) modified 
by the addition of the two second-order non-adiabatic 
potentials W + V.^5. . Further arguments in Ref. (22) are 
used to justify a certain normalization of the scattering 
function, but in essence one next introduces a new function 
X by the definition 


(1- r N ) 


(34) 
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and. a differential equation for X is derived "by direct 
substitution. Again neglecting terms higher than second 
order, one obtains the equation 


- [V, > G > 


( 55 ) 


where V is the "distortion potential" <(V G) 2 > appearing 

O X 

in Appendix A. From Eqs. (Ah) and ( 10 ) we find that _ 1 
at x - 0 . Thus, as Callaway et al 22 have emphasized, the 
distortion potential cancels the polarization potential V 2 
at the origin, and hence its inclusion seems to serve the 
same purpose as does monopole suppression 1 . Unfortunately 
Vp is fairly long range, decreasing like x 6 for large x } 
while V^q exponentially decreases. The equivalence of these 
two ways of suppressing excessive attraction is limited to 
small x, and consequently Eq. (35) gives far too much re- 
pulsion; the scattering length is^ 2 a - - O.783, even worse 
than the result of Eq. (32). 

k . We have previously reported loa a method of obtaining lower 
bounds to the phase shifts, which is simpler than that of the 
present work and not much different in its results. In 
essence it is a hybrid between the present method, which 
independently varies two functions and obtains coupled dif- 
ferential equations, and the method of Eq. (51) which varies 
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only one function. 

In this method, one again uses a trial function of 
the type represented by Eqs . (l) and (4), with 


F04)-= 


(36) 


where T(x) is some function which approaches unity for 
large x and decreases for small x. It is to contain 
adjustable constants 23 , and is designed to suppress the 
short-range correlation. Variation on X leads to the 
following equation, analogous to Eq. (31) : 


/ i . h - e! [-1 7 ; = o . 


(37) 


The optical potential of Eq. (32) is modified, and the 
lower bound phase shifts are solutions of the differential 


equation 
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v. 



*T*(V a **0 * (t 2 V n -2-VT'W)^ 


+ (JTV Z \ *) N 1 



( 38 ) 


which is seen to reduce to Eq_. ( 32 ) if T _ 1. For k - 0, 

L = 0 scattering, the optimum form of ’T would he that 
shown in Fig. 1, where T _ F/x, while F and X were numerically 
obtained. Fork^ 0, however, X and F both are oscillatory but 
their zeroes do not generally coincide, and one must use a 
smoothed approximate form for T, as has been done for the 
annihilation problem in Section III. It is also possible to 
combine the monopole suppression method with the present method. 
In that case, the constant a would measure the amount of mono- 
pole correlation retained, while the function T would be capable 
of making adjustments in the remaining parts of the correlation 
function . 

5. Two interesting modifications of the close-coupling approximation 
have been introduced recently. They are quite close in 
approach to the present coupled-equation method. 

a. Damburg and Karule 24 , noticing that the close -coupling 
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method fails to give the correct asymptotic dipole 
polarizability, have proposed a trial function which 
is equivalent to our Eqs. (l) and (4), with a dif- 
ferent form for G: 


C- tK (4 4) - (/i- + r A 0 fi 0**- *0 , 

(^CCr&- K. - A.® * ) 


(39) 


The motivation for this choice is clear: comparing 

Eq. (39) with Eq. (A8), one sees that G^ correctly 
describes the asymptotic form of the adiabatic cor- 
relation function, as far as the coordinates of the 
electron (r,T)) are concerned, and since F(x) can be 
freely varied, the solutions of the resulting coupled 
equations will be identical, for large x, to those 
of Eq. (8). Eq. (39) is separable in r, T)* x, however, 
and hence does not give the same short-range dipole 
correlation as does as does the Dalgarno-Lynn function G. 
In addition, each multipole term in the correlation 
function would require another coupled equation (as is 
the case for ordinary close -coupling) . In return one 
gains enormously in the simplicity of the potentials which 
must be used; all the potentials in this method can be 
evaluated analytically, while some reduce trivially to 
constants. The method is promising, but no numerical 
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results are available yet. 

b . Perkins 25 has made a modification starting from 
a conventional close-coupling expansion containing 
Is - 2p - 3d. states. Noting that, since only one 
term of each angular symmetry is retained^the ortho- 
gonality of closed and open channels [Eq. (2a)] is 
assured, he has modified the forms of the 2p and 3d 
functions as follows 


-^e a Pl , 

<k' A* e P* • 


0K>) 


The ususal coupled differential equations of close - 
coupling are solved numerically, but now two additional 
non-linear parameters, § and T), are adjusted to maximize 
the phase shifts. The results are much better than 
the corresponding unmodified close -coupling results, 
although the exact asymptotic polarizability is not 
achievable. The method can be extended to include higher 
i-values with the usual difficulty that more coupled 
equations are then required. Without generalizing the 
basis functions considerably, one would not expect the 
procedure to converge toward the exact wave function, but 
the lower bound theorem holds and offers the usual estimate 


of quality. 
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6. An interesting general formation was recently presented 26 , 
which uses numerical solutions of the two-center problem 
in the adiabatic approximation, and which may prove 
applicable to positron scattering. Similarities with 
our method include generalizations of our and W. The 
binding energy of the e + - e - e + molecule was computed 
as an example, but no scattering results were given. 
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Appendix A. Evaluation of the Potentials 

The analytic solution of Eq_. ( 5 ) was given by Dalgarno and Lynn 5 
and has the form 


G- = 

+• (_i +■ G- C x ) , 


»here f\(p£) -=. J |_ (*-* 0 ^ ~ 1 1 , 

G. (x-sjo) - El -JU (.\-s) 

- e’ 25t [^ (A l^ +< l) -^n. 0 yd]^ 


CD 


and E^(-Z) 


dy e ^/y. In terms of these elliptical coordinates, 
Z 


r = g x 
x and r. 
form 


+ p) • and cos 7] = (l + \p)/(\ 

The potential V _ — [l + — ] 

- x * p - \ 


+ p), where T| is the angle between 
? the bracket notation takes the 


oo 1 




X- 


l ix S e xtXv) (xV) -P (*,y) 


(A2) 
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and we may determine G(x) by requiring <G> _ 0, obtaining 

&(*•) ^)]+ e ^ X l Z(l Tk)XnY " («) 

where £ny = 0,5772 , Euler's constant. The solubility of Eq. ( 5 ) 

depended on its additive separability in elliptical coordinates. For 
the same reason, it is possible to evaluate most of the potentials 
defined in Eq. (9) in terms of one -dimensional integrals. Since, how- 
ever, the important function W cannot be reduced in this way, we have 
numerically evaluated all the potentials uniformly by the use of 
Gaussian quadrature in two dimensions . We have confirmed the accuracy 
of this procedure in three ways : We have evaluated <V>, <G>, and <GV> 

in elliptical coordinates, using our numerical method, and compared the 
results with the known values, obtaining excellent agreement. In particular, 
<G> never differs from zero by more than about 10 7 , and we feel that this 
figure measures the accuracy of the two-dimensional integrals. 

The two non-adiabatic terms, and W, present additional problems. 

To evaluate W directly would involve the operation V 2 G, which would 

X 

give a very intricate analytic expression. Instead, we have re-written 


the expression for W in the form 
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V = - z 








(A4) 


where V n = <(V G)^ is the "distortion potential" employed in Ref. 22, 
and discussed here in Section V. The function involves both x and 
n components of the gradient (with r held fixed), while only requires 
the x component. Both components were found analytically and inserted into 
the numerical double integration, while the first derivative of needed 
in Eq. (a 4) was gotten by numerical differentiation of V . All three 
terms in Eq. (kk) have asymptotic inverse-power behavior, with x” 6 as 
the leading term; these all cancel leaving short-range (exponential) 
terms dominant. 

In Section III we require the value of G(x, x), ie, G(r, x) evaluated 
at r _ x. In elliptical coordinates this relation becomes \ _ p = 1, 
since then r _ x and cos 1] - 1. Care is required in evaluating Q(x, l), 
since its first two terms are singular. Letting p = 1 - e and expanding 
the exponential integral about e _ 0, we have 
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C x (*, ' -4) « U (re x) - -JU e - e “ C 2X ) -i**] 


(A5) 


and the limit e - 0 can be taken. The final result is 


&(*,*)*= (uktlU*™-) + e 


1 - I Ei (-2^) +e 2aL ( J ^ x ) 


(A6) 


This is the function which, properly scaled, is used in the e + - He 
annihilation and forms the basis for Table III. 

To obtain the long-range asymptotic forms shown in Eq. (ll), we 
expand V in Legendre polynomials (for x > r) : 



-k I (xf 


(AT) 
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Then a particular solution of Eq. (5) is (for x > r) 



(A8) 


which is the correct asymptotic form of G. By inserting these expressions 
into the definitions, Eq. ( 9 ), one can carry out the integrations in 
spherical coordinates and obtain the leading terms, (in x 1 ) for each 
potential. Since V has no long-range & - 0 term, V x is di ort-ranged. 

? - ( 0 i 1 ) f \ 

To show that W is also, one notes that V x + P (cos 71) - 0, 
for 1 > 0. To illustrate, let us compute the asymptotic form of N: 



or 


|S/ 4l/(sx 4 ). 


(A10) 
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*» 


The expansion can be carried to any desired order in i, and represents 
an asymptotic expansion, since we have set the upper limit of the integral 
equal to oo rather than x. Care is required in evaluating the asymptotic 
form of V 3 since more than one value of a contributes to the same order 
in x" 1 , because three Legendre polynomials appear. 

The short-range behavior of the potentials (Eq. 10) is due to the 
small - I parts of the distortion function G. In Ref. 1. the i - 0 solution 
of Eq. (5) is exhibited. Near x = 0 it becomes 

G- 0 *Z J- -A' +- ! -X. a (_1-3L) (zr/l) (A11) 


and 


-2«) [A / zw ) M. - ~ 



(A12) 


Since, for small x, all the integrals have limits of 0 and oo, one finds 
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K> = <<**>« 

K - “ <& ^ G- 0 > 




Hoc. . 


(A15) 


The constant terms in Eq. (A13) are unaffected by the neglect of 
higher ^-values, but the second term in N and, hence, the leading 
term in V^, is modified by the l-l part of the expansion. One finds 


G, ~ Pi \v- * C^)J ^ 


(AHO 


1 2 

and, N 1 = <G^> — — x . No higher terms contribute to order x 1 2 , and 
we obtain the results given in Eq. (lO). The i _ 0 and l _ 1 parts of 


can now be found from Eq. (Ak ) . The result is ~ 0 and V 
near x = 0, and no higher terms contribute. 


D1 
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Appendix B. Solution of the Scattering Equations 

The numerical integration of Eq. (8), in partial wave decomposition, 
would be straightforward were it not for difficulties connected with the 
asymptotic form of the closed -channel part of the trial function, Eq. (l). 
The problem is very similar to those arising in the close -coupling method 
with closed channels, and our treatment follows the review by Burke and 
Smith 27 fairly closely. 

Since we are dealing here with a set of two coupled second-order 
differential equations, four linearly independent solutions can be found, 
with each solution consisting of a pair of functions (x, F). We will 
see later that only two of these solutions are sufficiently regular at 
x - 0 to be admissable. The general solution of Eq. (8) is thus a 
linear combination of these two solutions, and our problem is to determine 
such a combination which is sufficiently regular as x -* <». Making partial 
wave expansions as follows 

<50 . _ 

F (*) - £ oc ^( x ) f^C CoA - 6 ’'', (B1 

JL-6 


one can write the partial wave equivalent of Eq. (8) : 
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(B2a) 


« 3 » 4 . 

- M' 1 V w - V ’ V, 4 . 


(B2b) 


From Eq. (10) we obtain the dominant terms of the potentials near x _ 0. 
Representing each linearly independent solution as a vector ^ = (U, g), 
we find two linearly independent regular solutions near zero: 


<fc*) = (^>°), 
( x ) - 


where t^ « x^ + 1 [1 + x/(L + l)], plus higher powers of x. Eq. (B2) 
is numerically integrated, using these starting forms, up to a point 
x _ x, and the general solution there is 


. ri>) 

ip L (s.\ - A f L (5.) +■ 8 (|l 


(B4) 
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where A, B, are constants to be determined. For large distance 
we can re-write Eq. (B2) using the asymptotic expression in Eq. (ll). 
Retaining terms which fall off less rapidly than x 3 we find 



(B5a) 


n 




‘i. a' 




A?> u ‘ 


(B5b) 


[We will discuss the L = 0, k _ 0 case below.] Taking x ^ 20 as a 

typical "asymptotic" distance, we see that the x 2 term is not neglibible 

for small k( !=a 0.1) in Eq. (B5a), but is negligible in Eq. (B5b) . Of the 

four solutions to Eq. (B5), three are well-behaved for x -» <». If we 

set U - 0 in Eq. (B5b) and neglect x 2 terms, we find that the closed- channel 

function g^ has two solutions, 

5± e e 1 *■ ^ 




(b6) 



The asymptotic regularity condition, Eq. (2b), requires 


o° 

(*)<*>, < B 7> 

J o ' J 

— A 3/2 

and since N ~ x , must increase slower than x for large x. 

Thus we are forced to drop the rising exponential form S+ of Eq. (B6) . 

The non-zero linearly independent solutions of Eq. (B5a) are composed 
of the spherical Bessel functions : 


U u a x , « 

= oc fii L . 


(B8) 


Corresponding to these open-channel solutions one finds a slowly 
converging 28 series of inverse powers for g^. The three solutions 
for large x are thus : 


(*) X (x m u) 


where 
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, oc 1 



The coefficients in Eq. (BIO) can be obtained, directly by use of 
Eq. (B5b) and equating like powers of x. The result for L - 0 is, 
for example. 







4 



(Bll) 

L 


Eq. (B2) is numerically integrated inward from some "asymptotic" point 
x = X o’ us ^ n S “these starting forms, down to the matching point x — x, 
where the general solution is 



- t3 - 



6b 
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Setting ^ Q - ‘H (&> and 



gives four linear equations in the unknowns A, B, C and tan 6^. 

We have tested the stability and consistency of the numerical procedure 
by comparing the results obtained for x q = 20, 25 , 30 and x - 2, 3> 4, 
and find essentially no variation in the phase shifts for L = 0 (see 
Table IV). (An empirical correction for L - 1, obtained from our 
previous non-variational work was applied to account for the increase 
produced by the long-range x 4 potential beyond x - 20. For very low 
energies exact analytic expressions exist 29 and are superior to our 
numerical results} 

For L = 0 and k = 0 the long-range potential V 2 is not negligible in 
Eq. (B2a), even for x 20, so it is not correct to use the asymptotic 
form Eq. (B5a) . Instead, one has 



(Blja) 



- 44 - 


3 


ii 

c 


4 Q t _ 3^_ 

5T 45 


(B15b) 


Solution is unchanged from the form given in Eqs. (B9) and (B6), 

but the other two asymptotic solutions are now 





36 oc 


j 





'foe- 


(B14) 


and the linear combination shown in Eq. (B12) is unchanged except that 
tan 6 Q is replaced by the scattering length a. 
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X 

V ! 

-V 2 

v 3 

N 

V N 

W 

0 

*a-o 

1 

IG- 1 ) 

3/4 

10 

T x 

6(1 - 2 x ) 

0.5 

2.207 

6. 96( -1) 

1.229 

5-38(-i) 

5-54(-i) 

1.913 

1.0 

5-4i(-i) 

3-99 

i-57(-i) 

3.13 

3.34 

6.42(-l) 

1-5 

1.66 

2.34 

-1.40(-2) 

1.89 

1.79 

2.62 

2.0 

5-49(-2) 

1.41 

-3-99 

1.20 

1.07 

1.17 

2-5 

1.89 

8.48(-2) 

-3.32 

7 - 73( -2) 

6.84(-2) 

5-25(-2) 

3.0 

6.6l( -3) 

5.14 

-2.18 

4.96 

4.41 

2.28 

3-5 

2.34 

3-14 

-1.29 

3.19 

2.80 

9- 61 (-3) 

4.0 

8.39(-4) 

1-95 

-7-20(-3) 

2.06 

1.77 

3-96 

4.5 

3.02 

1.25 

-3.91 

1.36 

1.11 

1.61 

5.0 

1.09 

8.19(-3) 

-2.11 

9-l2(-3) 

7 . 04 ( - 3 ) 

6.43(-4) 

6.0 

1.43(-5) 

3.86 

-6. 23 (-4) 

4.43 

2.98 

1.00 

7.0 

l.90(-6) 

2.03 

- 2.00 

2.36 

1.38 

i-53(-5) 

8.0 

2-53(-7) 

1.17 

-7-l8(-5) 

1-37 

7.00( -4) 

2.3l(-6) 

9.0 

3 . 38 ( -8 ) 

7-l8(-4) 

-2.91 

8.47( -4) 

3.83 

3-43(-7) 

10.0 

4.54(-9) 

4.67 

-1.31 

5.52 

2.24 

5.06(-8) 

TABLE 

I. Numerical 

values for the potentials defined 

in Eq. (9), 

evaluated 


using the 

complete 

form of G in 

elliptical 

coordinates . 

(Negative 


exponents 

appear in 

parenthesis 

in the usual way) . 
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E(ev) 

A 

B 

0 

I.83 

3.16 

0.54 

I.65 

2-79 

1.22 

1-57 

2.58 

2.17 

1.51 

2.43 

3.40 

1.49 

2-35 

4.90 

1.49 

2.31 

6.66 

1.52 

2.30 

8.70 

1-57 

2-33 

11.0 

1.62 

2.38 

13.6 

1.68 

2.44 

17.8 

1.78 

2.54 


TABLE II. Positron -Helium Annihilation Rate — . Column A 

gives the present results, and 001111121 B is the 
non-variational result of Ref. 15 - The Dirac rate 

is K = 
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V = £x 

M 

Y 

0 

1.67 

1 

0.5 

1.68 

1.91 

1.0 

1.83 

2.78 

1-5 

2.13 

3.50 

2.0 

2.55 

4.07 

2-5 

3-07 

4.53 

3.0 

3 • 66 

4.92 

3-5 

4.29 

5.25 

4.0 

4.90 

5-55 

4.5 

5-47 

5.82 

5.0 

5-97 

6.06 


TABLE III. Correlation functions for e + - Helium annihilation. 

M is used in the present work and is defined in Eq. (28). 
Y = [l + G(x, x ,)] 2 + <G 2 > is used analogously in the 
non-variational work (Ref. 15) , and employs full monopole 


suppression in G. 
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tan § (k k O.l) 


X 

o 

OJ 

II 

o 

X 

IfN 

OJ 

II 

X° 

X 

O 

II 

V-M 

O 

2 

0.12753 

0.12804 

0.12808 

5 

0.12747 

0.12798 

0.12803 

4 

0.12675 

0.12726 

0.12731 


X 

0 

OJ 

II 

0 

X 

tan 6 (k = 0 - 7 ) 

x _ 25 
0 

x - 30 
0 ~ 

2 

- 0.07393 

-O.O7386 

- 0.07383 

3 

-O.O7396 

-0 . 07389 

-O.07386 

4 

- 0.07440 

- 0.07433 

- 0.07430 


TABLE IV. S-wave e + - H results. The sensitivity of tan 6 to 

matching radius (x) and asymptotic radius (x q ) is shown, 
for two values of positron momentum k. 
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k 

tan 6 

tan fi 

s 

tan 6 h 

Q 

Sis 

0 

-1.85 

-2.10 

O.582 

0.91 


0.1 

0.128 

0.152 

-0.058 

0.89 


0.2 

O.I58 

0.190 

-0.116 

0.90 

Q.87 

0-5 

0.135 

0.170 

-0.170 

0.90 


0.4 

O.O89 

0.121 

-0 . 222 

0.91 

O.89 

0.5 

0.034 

0.062 

-0.270 

0.92 


0.6 

-0.022 

0.007 

-0.314 

0.91 

O.89 

0.7 

-0.074 

-0.054 

-0.354 

0.93 



TABLE V. S-wave e + - H scattering results. The entries without 

subscripts refer to the present work, and are averages over x. 
The other phase shifts are: 6 —Schwartz (Ref. 6); Hartree 

S xi 

and 6gg— Hahn and Spruch (Ref. 16) . The figure of merit Q is 
defined in the text. For k _ 0, scattering lengths are given. 


- 53 - 


k 

tan 5 

tan 6 h 

fan 6 W 

Q 

0 

-0.511 

0.420 

-0 . 659 

0.86 

0.1 

O.O56 

-0.042 

0.050 

0.85 

0.2 

0.047 

-O.O85 

0.072 

0.84 

0.3 

O.O59 

-0.124 

0.071 

0.84 

0.4 

0.020 

-O.I65 

O.O56 

0.84 

0-5 

-0.007 

-0.205 

0.052 

0.85 

0.6 

-0.059 

-0.240 

0.002 

0.85 

0.7 

-0.075 

-0.278 

-0.051 

0.85 

0.8 

-0.107 

-0.513 

-0.066 

0.85 

0.9 

-0.142 

-0.546 

-0.100 

0.85 

1.0 

-0.176 

-0.579 

-0.154 

0.85 

1.1 

-0.208 

-0.409 

-0.168 

0.85 

1.145 

-0.225 

-0.425 

-0.182 

0.85 


TABLE VI. S -wave e + - He scattering results. The entries without 

subscripts refer to the present work and are averages over x. 
The other phase shifts are: 6 H - Hartree ; 5^ - nonvariational 
(Ref. 2). For k = 0, scattering lengths are given. Q is 
computed assuming 6^ to be exact. 
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k 

tan 

tan 6 X (LB) 

tan 6 -l (ext) 

0.1 

0 . 0073 ( 1 ) 

0.0082 

0 . 0086 ( 1 ) 

0.2 

0 . 0263 ( 2 ) 

0.0289 

0 . 032 ( 1 ) 

0.3 

0 . 0518 ( 5 ) 

0.0548 

0.066(4) 

0.4 

0.0764(9) 

0.0801 

0.11(1) 

0.5 

0.0952(14) 

0.0994 

0.14(1) 

0.6 

0 . 106 ( 2 ) 

0.112 

0 . 17 ( 2 ) 

0.7 

0 . 109 ( 2 ) 

0.119 

0 . 19 ( 2 ) 


TABLE VII. P-wave e + - H scattering results. The present results 

(corrected for long-range effects as described in 

Appendix B) are compared with the lower hounds (LB) 

and extrapolated Values (ext), obtained in Ref. 21. 

The estimated error in the last digit is shown in 


parenthesis . 
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Figure Caption 

Figure 1. Zero-energy scattering solutions for e + - H. The 

functions X and F are open -channel and closed -channel 
functions, and T _ F/x is the suppression factor. The 
asymptotic normalization is X = 1 - a/x, F _ X and 


T _ 1. 
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